Generalised parton distributions at HERA and prospects for COMPASS 
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We present a model of generalised parton distributions based on a forward ansatz in the DGLAP 
region. We discuss some aspects of the parametrisations, as the dependence in t, with factorised and 
non-factorised approaches, where t is the square of the four-momentum exchanged at the hadron 
vertex. The predictions of this model are then compared with DVCS cross sections from HI and 
ZEUS, and a related observable, the skewing factor, defined as the following ratio imaginary ampli- 
tudes : R = Im A (7* + p — » 7 + p)\ /Im A (7* + p — » 7* + p)\ . It is an interesting quantity 

lt=o lt=o 
including both the non-forward kinematics and the non-diagonal effects. Finally, we discuss the 

beam charge asymmetry, which is certainly the most sensitive observable to the different hypothesis 
needed in the GPDs parametrisations. We show that the approximations done for the t dependence 
lead to significant differences for the predictions in the HERMES kinematic domain and prospects 
are given for COMPASS. 
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I. INTRODUCTION 



Measurements of the deep-inelastic scattering (DIS) of leptons and nucleons, e+p — » e + X, allow the extraction of 
Parton Distribution Functions (PDFs) which describe the longitudinal momentum carried by the quarks, anti-quarks 
and gluons that make up the fast-moving nucleons. These functions have been measured over a wide kinematic range 
in the Bjorken scaling variable xbj and the photon virtuality Q 2 . While PDFs provide crucial input to perturbative 
Quantum Chromodynamic (QCD) calculations of processes involving hadrons, they do not provide a complete picture 
1— 1 of the partonic structure of nucleons. In particular, PDFs contain neither information on the correlations between 
i-C ' partons nor on their transverse motion, then a vital knowledge about the three dimensional structure of the nucleon is 
lost. Hard exclusive processes, in which the nucleon remains intact, have emerged in recent years as prime candidates 
O ■ to complement this essentially one dimentional picture 0, i, i, H B B 0, H S El, El El El, El- This missing 
information is then encoded in Generalised Parton Distributions (GPDs). These functions carry information on both 
the longitudinal and the transverse distribution of partons. The recent strong interest in GPDs was stimulated by 
their relation with the spin structure of the nucleon. Indeed, GPDs are so far the only known means of probing 
[ — . the orbital motion of partons in the nucleon through Ji's Sum Rule [l[ , which relates unpolarised GPDs to the total 
^ angular momentum of the proton. 

The simplest process sensitive to GPDs is deeply virtual Compton scattering (DVCS) or exclusive production of 
real photon, e+p — > e + j+p. This process is of particular interest as it has both a clear experimental signature and is 
calculable in perturbative QCD. Also, it does not suffer from the uncertainties caused by the lack of understanding of 
the meson wave function that plague exclusive vector meson electroproduction. The DVCS reaction can be regarded 
as the elastic scattering of the virtual photon off the proton via a colourless exchange, producing a real photon in 
the final state. In the Bjorken scaling regime, QCD calculations assume that the exchange involves two partons, 
having different longitudinal and transverse momenta, in a colourless configuration. These unequal momenta are a 
consequence of the mass difference between the incoming virtual photon and the outgoing real photon. The DVCS 
cross section depends, therefore, on GPDs 0, !, i, @, i, M, El El, El, 0] • 

The reaction e+p — > e + 7 + p receives contributions from both the DVCS process, whose origin lies in the strong 
interaction, and the purely electromagnetic Bethe-Heitler (BH) process where the photon is emitted from the positron. 
Therefore, measurements can be done either of the DVCS cross section itself or of the interference between the DVCS 
and BH processes, accessible via asymmetries. Measurements of the DVCS cross section at high energy have been 
obtained by HI [3, El] and ZEUS [l7| and the helicity asymmetry in DVCS has been measured at lower energy with 
polarised lepton beams (and unpolarised targets) by HERMES [l8( and CLAS [19|, l2fj . Recently, lepton beam charge 
asymmetry have been measured by HERMES 21]. These last observables are very sensitive to the shape of the GPDs 
and give a good opportunity to pin-down their parametrisations. It is the goal of a future program at COMPASS 
pi. I22I] to realise a dedicated measurement and to access the beam charge asymmetry under favorable conditions. 
In this paper, we focus our study on the HERA results and we show that GPDs parametrisations based on the 
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model of Ref. [Lj, [lj] are well suited to describe all measurements in this range of xsj {^Bj < 0.1). In section II, we 
present the model. We show the predicted DVCS cross section compared to data in section III. In this section, we 
also extract the skeewing factor from data and from the model. This factor is an interesting observable to illustrate 
the data/model comparison and the building of the skewedness. In section IV, we discuss the case of the beam charge 
asymmetry measured at HERMES and we give some prospects for COMPASS. 



II. MODEL OF GENERALISED PARTON DISTRIBUTIONS 



As mentioned in the introduction, GPDs are an extension of the well-known parton distribution functions (PDFs) 
appearing in inclusive processes such as deep inelastic scattering (DIS), or Drell-Yan, and encode additional infor- 
mation about the partonic structure of hadrons, above and beyond that of conventional PDFs. As such, the GPDs 
depend on four variables (X, £, Q 2 ,t) rather than just two (xsj, Q 2 ) as is the case for PDFs. We use the notations of 
Ref. [2|, where the variable X is defined in the range [0, 1] and the variable £ defines the skewedness with £ ~ XBj- 
This allows an extended mapping of the dynamical behavior of the nuclcon in the two extra variables, the skewedness 
£, and the four-momentum squared exchanged at the hadron vertex, t. In the following, we define GPDs at a starting 
scale Qq and their Q 2 evolution is generated by perturbative QCD [3, HH [HI, EH 0, [3] • They contain, in addition to 
the usual PDF-type information residing in the so-called "DGLAP region" (for which the momentum fraction variable 
is larger than the skewedness parameter, X > £), supplementary information about the distribution amplitudes of 
virtual "meson-like" states in the nucleon in the so-called "ERBL region" (X < £) [9] . The parametrisation of GPDs 
is then a complicated task as the degrees of freedom are much larger than for PDFs. 

Concerning the A, £ dependences, we follow the model presented in Ref. [H, [Hj], using a forward ansatz at the 
initial scale Qq = 1.3 GeV for the DGLAP domain [£, 1] : 



X-C/2 
l-C/2 



g *(*0 g l-C/2 ' (1) 

where qs refers to a singlet forward distribution [23[ and Hs(X, £) = H q (X, £) + Hq(X, £). A similar relation holds for 
valence and gluon distributions. This Ansatz in the DGLAP region corresponds to a double distribution model with 
an extremal profile function allowing no additional skewedness except for the kinematical one. Concerning the ERBL 
domain, we follow the strategy of Ref. [HI [HI • GPDs must be continuous at X = £ and verify correct symmetry 
relations at X = £/2. Then, the singlet distribution can be parametrised as 

H S (X, £) = ffs(C, 0^/^[l + As(0f(X, £)] (2) 

where f(X, £) is an even function of the variable X — £/2 such that /(£,£) = to ensure the continuity relation at 
X = £. Thus, we can write : /(A, £) cx 1 — (X — £/2) 2 / (£/2) 2 . The parameters As(£) are determined for each value 
of £ by requiring that the first momentum sum rule of the GPD is satisfied [2, 0] • All other GPDs are parametrised 
as exposed in Ref. [HI, EH- The skewed QCD evolution code has been rewritten for the purpose of this paper. 

Results are shown in Fig. [T] for the generalised u quark singlet density, for two values skewedness £ — 0.05 and 
£ = 10 -3 , at the initial scale and after a QCD evolution. We notice that the skewedness effect manifest itself 
much before the border between the DGLAP and ERBL domains at X = £. In particular, for the lowest £ value 
of 10 -3 , the GPD follows exactly the PDF till X ~ 10~ 2 , where the skewing effect starts to appear, as expected 
from the forward limit condition : GPDs are obliged to reproduce the PDFs in the forward limit £ — + 0, with 
Hjj s (X, £ = 0) = Us(X) = u(X) + u(X). Also, we observe in Fig. [T]that the generalised u quark singlet distribution 
is antisymmetric at X = £/2 in the ERBL region and this property is preserved under evolution. In Fig. [21 we present 
the generalised gluon density again for two values skewedness ( = 0.05 and £ = 10 -3 , at the initial scale and after a 
QCD evolution. Here also, for the lowest £ value of 10 -3 , the GPD follows exactly the PDF till X ~ 10~ 2 as expected 
from the forward limit condition Hq{X, £ = 0) = XG(X). 

Concerning the t dependence in the low XBj kinematic domain of HI and ZEUS measurements (l6j . it has been 
shown that the DVCS cross section, da/dt, can be factorised and approximated by an exponential form e _fl '*' , implying 
a factorised dependence in e~ h l 2 for GPDs. It follows from Ref. [n| that the measured i-slopes b present a Q 2 
dependence compatible with the form : b(Q 2 ) = A (l — 0.15 log(Q 2 /2)) , where A = 7.6 ± 0.80 GeV -2 . We keep this 
expression in the following when describing HI and ZEUS measurements, as it is done also in Ref. JTl|. Of course, it 
can only be considered as an effective approximation of a more realistic non- factorised approach [fl^In Ref. [l]| . a 
non-factorised t dependence is proposed with a Regge motivated approach. In this approximation, we can write the 
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FIG. 1: The generalised distribution of the u quark singlet (multiplied by X) XH Ua (X,() = XH U (X, () + XH U (XX) is 
presented for two values of the skewedness £ = 0.05 and 10 -3 . The standard PDF XUs(X) = Xu(X) + Xu(X) is also 
displayed. These functions are given at the initial Q% value of 1.69 GeV 2 (left) and after the NLO QCD evolution towards 7.14 
GeV 2 (right). In case of GPD, the evolution is done in the DGLAP and ERBL domain ensuring proper continuity relations 
at the border [H, Two vertical lines at X = 0.05 and X — 10~ 3 , corresponding to the skewing values, are also shown to 
indicate the position in X of the transition between the DGLAP and ERBL domains. 



singlet distribution at the initial scale in the DGLAP domain : 

l-C/2 



H s (X,(,t;QZ) 



X-(/2 



qs 



X-C/2 
l-C/2 



l-C/2 



(3) 



with similar expressions for the gluon and valence distributions. The ERBL domain is parametrised following Eq. ([2|) 

l-C/2 Q ' s * 

to take into account the t dependence. The continuity at the border X = £ is then 



with the global factor x-c/2 

ensured. This parametrisation is mixing the X and t dependences, thus it is labeled as non-factorised. It is insprired 
from Regge phenomenology as parameter a' s is characteristic of a Regge trajectory. We take the values : a' s = 0.9 
GeV" 2 and a' g = 0.5 GeV -2 as in Ref. [HI- 

In the next sections, we compare this model with the different approximations on the t dependence to available 
observables provided by HI, ZEUS and HERMES experiments. 



III. PREDICTIONS FOR THE DVCS CROSS SECTION AND RELATED OBSERVABLE 

Using these paramctrisations of GPDs at the initial scale Qq = 1.3 GeV, calculated to higher Q 2 values with a 
skewed QCD evolution, and assumming the t dependences described above, we can provi de p redictions for the DVCS 
cross section measurements. We present the results in Fig. [3] with HI and ZEUS data jl6l. fl7l|. The measurements 
are integrated over t in the range [—1,0] GeV 2 . Then, the model predictions are also integrated over the same range. 
Concerning these predictions, the uncertainty due to the limited knowledge of the ^-slopes is displayed in Fig. [31 The 
calculations in the non-factorised t approach fall within this error, which shows that the DVCS cross sections are not 
sufficient to discriminate between these two hypothesis. The same conclusion holds if we consider the DVCS cross 
section differential in t, dcr/dt, measured in Ref. [la ]. Indeed, the DVCS cross sections are senstive essentially to the 
values of the GPDs at X = ( ~ xsj (at low xsj < 0.01 for HI and ZEUS kinematic range). Thus, the global factor 




FIG. 2: The generalised distribution of the gluon Ha(X,C) i s presented for two values of the skewedness £ = 0.05 and 10 3 . 
The standard PDF XG(X) is also displayed. These functions are given at the initial Qq value of 1.69 GeV 2 (left) and after 
the NLO QCD evolution towards 7.14 GeV 2 (right). 



a' t 

of the non factorised approach ^ x - c/2 ) can ^ e wl "ittcn as cxp(— a' ln(10~ 3 )|i|) for HI kinematics [16], which is 

reproducing to a very good approximation the factorised dependence in exp(— b\t\) with b ~ 6 GeV -2 . 

Also, a good agreement between data and the model predictions is observed : it shows that the simple description 
of the initial condition using a forward ansatz in the DGLAP domain and then generating the higher Q 2 values from 
a skewed QCD evolution is a good approach. In particular, it can be noticed in Fig. [3] that the Q 2 dependence of 
the DVCS cross section is well reproduced. The W dependence can be fitted by a form in W s with 8 ~ 0.7, a large 
exponent value characteristic of a hard QCD process [m. \vf\. 

It is interesting to present the DVCS cross section measurements in a form which illustrates directly the skewing 
properties by defining a quantity which includes both the non-forward kinematics and the non-diagonal effects. 
Namely, we set the ratio between the imaginary parts of the DIS and DVCS (forward) scattering amplitudes at zero 
momentum transfer 1241: 



Im A (7* + p — > 7 + p) 

R = 



t=o 



Im A (7* + p — ► 7* + p) 



(4) 



t=o 



The virtual photon is assumed to be mainly transversly polarised in the case of DVCS process due to the real photon 
in the final state and has to be taken transversly polarised in the DIS amplitude also. The scattering amplitude for the 
DIS process can be obtained from the (transverse) DIS cross section [25|, that is ImA{^f*p — > 7*p) ~ <jt{"1*P — * X). 
In fact, the DIS amplitude can be written using the usual pQCD fits for the transverse proton structure function, 
ctJ p — (4ir 2 aEM IQ 2 ) Ft( x 'Q 2 ) anc ^ * ne DVCS scattering amplitude can be obtained from the recent measurements 
on the Deeply Virtual Compton Scattering cross section. We have shown that the DVCS cross sections are not 
sensitive to the hypothesis done on the t dependence described in section II. Then, obviously, the same conclusion 
holds for the skeewing factor and we keep, for its determination, the effective approximation supported by the HI 
data p3 that the t dependence of the amplitude can be factorised out and parameterised as an exponential. The 
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FIG. 3: DVCS cross section as a function of Q 2 at W = 82 GeV (left), and as a function of W at Q 2 = 8 GeV 2 (right). The 
GPD model predictions [13, [3 are integrated over the t range of the measurements and displayed (see text). The model of 
section II, using an exponential factorised t dependence is shown with full lines on both the Q 2 and W dependences of the 
DVCS cross section. The dotted lines illustrate the uncertainty on the t-slope dependence as described in section II. On the 
DVCS cross section as a function of W, we also display, as a dashed line, the model described in section II using a non-factorised 
t dependence (see text). The results of the factorised and non-factorised approaches are close and the differences are within 
the dotted lines. For clarity, we do not show the non-factorised curve on the Q 2 dependence plot, as it would be very close to 
the factorised model on a logarithmic scale. 



total DVCS cross section can be related to the corresponding amplitude at t = 

1 2 



ImA (7*p — * 7p) 



t=o 



16tt6(Q 2 ) 



(5) 



where b(Q 2 ) is the i-slope function of Q 2 described in section II. The expression above is corrected by taking into 
account the contribution from the real part of the amplitude by multiplying Eq. by a factor (1 + p 2 ), where p is 
the ratio of the real to imaginary parts of the DVCS amplitude. 

Following Ref. [2(|, a discussion concerning the calculation of the real part in the forward model is necessary. At 
the initial scale Qo = 1.3 GeV, the input GPD parametrisations in this model do not fulfill all polynomiality relations 
[U d(| as only the first momentum sum rule of the GPD is imposed (see section II and Ref. [ill, HH 1 )- In Ref- [III, 
a discussion of a simplified version of the forward model at the initial scale Qo is presented, which shows that the 
calculation of the real part of the DVCS amplitude is not correct at this scale as it does not verify the dispersion 
relations. The mismatch is arising mainly from a lack of flexibility of the skewedness dependence at this very low scale. 
However, as soon as we let a sufficient range in Q 2 for the QCD evolution to take place, the skewedness dependence is 
determined mainly from the QCD evolution itself [HI, EH- ^ n s i m ple words, it means that the memory of the original 
input dependence in the skewing variable is vanishing with the QCD evolution in Q 2 (27[. In the range of HI and 
ZEUS data [H,[]3, f° r Q 2 > 4 GeV 2 , we have checked that the real part amplitude is compatible with the dispersion 



1 A general procedure is also described in Ref. [l3t "T3 | to correct the input distributions in order to verify the polynomiality relations at 
all orders. In the future a new version of the model will be proposed applying these ideas. 



relations to a good approximation, better that 10 %. Then, we do not face the problem mentioned in Ref. [26j |. 
Namely, the ratio of real to the imaginary parts of the DVCS amplitude, p, can be calculated from the amplitudes 
determined in the model or using dispersion relations. In this last case, we can write p ~ tan(-7rA/2), where A = A(Q 2 ) 
is the effective power of the Bjorkcn XBj dependence of the imaginary part of the amplitude. Hence, in the range of 
HI and ZEUS data, for Q 2 > 4 GeV 2 , we use this property to correct the skewing factor extracted from the data. In 
this case, we estimate of the real part contribution by using the effective power for the inclusive deep inelastic reaction 
taken from Ref. (28J. For the theory prediction, we use the real part of the DVCS amplitude as derived following the 
model of section II. Finally, for the kinematic window available at HERA, the typical contribution of the term in p 2 
to Eq. © is of the order of 10 %. 

Considering the calculation discussed above, we can rewrite the skewing factor as a function of the cross sections 
for DIS (or) and DVCS : 



R = 4 y/ir apvcs b(Q 2 ) = y/<Jpvcs Q A KQ 2 ) 

* T (TP^X)V(1 + P 2 ) V^a EM F T (x,Q 2 )^/(l+p 2 ) ' 

Main theoretical uncertainties come from the t-slope, b{Q 2 ), given in section II. Results are shown in Fig. |4j where a 
good data/model agreement is observed within errors. 

On the skewing factor, we can exemplify the part of the skewing arising from the kinematic of the DVCS process 
and from the Q 2 evolution itself. Then, we apply the forward ansatz, used at the initial scale 1.3 GeV in the model 

qs ( x ~**/ 2 -Q 2 ) 

described in section II, at all values of Q 2 . It means that we impose the parametrisation Hs(X, (;Q 2 ) = — \~-£/2 — 
in the DGLAP domain for all values of Q 2 , and similar relations for valence and gluon distributions. As illustrated 
in figure HJ the measurements show that such an approximation, which only takes into account the kinematical 
skewedness, is not sufficient to reproduce the total skewing effects generated by the QCD evolution equations. 
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FIG. 4: Skewing factor R = Im A (7* + p — > 7 + p)\ /Im A (7* + p — » 7* + p)\ extracted from DVCS and DIS cross 

lt=o lt=o 
sections as explained in section III. The GPD model is also displayed and gives a good agreement of the data (full line). The 

forward ansatz model, used at all values of Q 2 , fails to reproduce the total skewing effects generated by the QCD evolution 

(dashed line) -see text-. 



IV. BEAM CHARGE ASYMMETRY (BCA) AT HERMES AND COMPASS 

The determination of a cross section asymmetry with respect to the beam charge. (da + — do~ ) / \da + + rfcr - ), has 
been realised by the HERMES experiment [U for x Bj ~ 0.1, Q 2 ~ 3 GeV 2 and \t\ < 0.7 GeV 2 . The interest of this 
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TABLE I: The cos <f> amplitude of the beam-charge asymmetry per kinematic bin in —t after background correction and the 
respective average kinematic values. In the last two columns, we give the predictions for the GPD model in the factorised and 
non-factorised dependences of t. 



measurement lies in its large sentivitythe shape of the GPD in X, £ [l3j, LLJ] and also in the correlations between 
longitudinal and transverse variables 

Using DVCS cross section measurements in the low XBj kinematic domain, we have tested with success the GPD 
model described in section II. In this section, we compare the predictions of this model with the results of HERMES 
[2lT |. which provides measurements of the cos(</)) amplitude of the BCA, where 4> is the angle between the plane 
containing the incoming and outgoing leptons and the plane defined by the virtual and real photon. In Table [H we 
compare the experimental values to the predictions of the GPD model described in this paper, in case of the factorised 
exponential t dependence and in the the Regge non-factorised t behaviour. A table is a good way to present these 
results as for each t values, the x and Q 2 values are different. As mentioned in the previous section, we need to 
let a sufficient range in Q 2 between the initial scale of the model (Qo = 1.3 GeV) and the measured values for the 
calculations of the real part of the DVCS amplitude to be completly correct. We have already discussed that for Q 2 > 4 
GeV 2 , the exact calculation of real part of the DVCS amplitude in the forward model and the derivation through 
the dispersion relation are well compatible. In the HERMES domain, for Q 2 ~ 2.5 GeV 2 , the agreement between 
both methods is of about 15 % and less good at lower Q 2 . Then, we consider only the HERMES measurements with 
Q 2 > 2.5 GeV 2 . 

Contrary to the DVCS cross section predictions, the calculations presented in Table U are very different in both 
approaches, factorised and non-factorised t dependences. Therefore, we confirm that BCA is a very sensitive and 
discriminating observable to study GPD models. Within the present large experimental errors, as shown in Table [H 
we can not favor one approach but it is clear that further high precision measurements at COMPASS would be of 
high significance (22J. 

In Fig. we compare predictions of our model to simulations of the BCA extraction at COMPASS using a muon 
beam of 100 GeV jJH2|- We present the comparison for one value of Q 2 (4 GeV 2 ) and two values of XBj (0.05 and 
0.1), as we know that for Q 2 > 4 GeV 2 , we stand in a safe kinematic domain concerning the real part calculations in 
the forward model (see above). When we compute the BCA in the factorised exponential t dependence approximation, 
we find values compatible with zero, which are not represented in Fig. 03 We just display the predictions of the model 
obtained in the non-factorised case. Then, both the cos(0) and cos(20) terms contribute to a significant level to the 
BCA at COMPASS, as illustrated in Fig. [3 We remark that our predictions do not match with the COMPASS 
simulation done with the model of Ref . |J| , which is another illustration of the large discriminative power of this 
observable on GPDs parametrisations. 



V. CONCLUSION 



An outstanding task in QCD is related to the extraction of GPDs. In contrast to PDFs, these functions contain 
informations on the correlations between partons, on their transverse motion and thus on the three dimensional 
structure of the nucleon. Exclusive production of real photon is a prime measurement to access the GPDs, either 
from DVCS cross section or from BH/DVCS interference. Data are getting more precise, leading to more refinement 
in the models. In this paper, we have addressed the case of the kinematic domain of HERA. We have shown some 
basic features, as skewing or Q 2 evolution of GPDs, in a model based on a forward ansatz in the DGLAP region. We 
have discussed the t dependence in the factorised and non-factorised approximations. A good agreement of DVCS 
cross sections from HI and ZEUS have been obtained in both cases. We have compared this model with success also 

on the skewing factor R = ImA (7* + p — > 7 + p) /ImA (7* + p — > 7* + p) , which is an interesting quantity 

t=o t=o 
as it includes both the non-forward kinematics and the non-diagonal effects. Finally, we have shown that the BCA is 

an interesting sensitive observable to the different hypothesis needed in the GPDs parametrisations. We have shown 
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FIG. 5: Simulation of the azimuthal angular distribution of the beam charge asymmetry measurable at COMPASS at = 100 
GeV. We present the projected values and error bars in the range \t\ < 0.6 GeV 2 for 2 values of XBj (0.05 and 0.1) at Q 2 — 4 
GeV 2 (see Ref. [II]). The prediction of the GPD model with a non-factorised t dependence is shown (full line) . The case of a 
factorised t dependence would lead to a prediction of the BCA compatible with zero and is not displayed. 

that the approximations done for the t dependence lead to significant differences for the predictions in the HERMES 
kinematic domain, and even larger for COMPASS. 
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